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ABSTRACT 



We constrain the form of the primordial power spectrum using Wilkinson Microwave 
Anisotropy Probe (WMAP) 3-year cosmic microwave background (CMB) data (+ other high 
resolution CMB experiments) in addition to complementary large-scale structure (LSS) data: 
2dF, SDSS, Ly-a forest and luminous red galaxy (LRG) data from the SDSS catalogue. We 
extend the work of the WMAP team to that of a fully Bayesian approach whereby we com- 
pute the comparative Bayesian evidence in addition to parameter estimates for a collection 
of seven models: (i) a scale invariant Harrison-Zel'dovich (H-Z) spectrum; (ii) a power-law; 
(iii) a running spectral index; (iv) a broken spectrum; (v) a power-law with an abrupt cutoff 
on large-scales; (vi) a reconstruction of the spectrum in eigh t bins in wavenumber; and (vii) a 
spectrum resulting from a cosmological model proposed by Lasenb v~& Doranl d2005l) (L-D). 
Using a basic dataset of WMAP3 + other CMB + 2dF + SDSS our analysis confirms that a 
scale-invariant spectrum is disfavoured by between 0.7 and 1 .7 units of log evidence (depend- 
ing on priors chosen) when compared with a power-law tilt. Moreover a running spectrum 
is now significantly preferred, but only when using the most constraining set of priors. The 
addition of Ly-a and LRG data independently both suggest much lower values of the running 
index than with basic dataset alone and interestingly the inclusion of Ly-a significantly dis- 
favours a running parameterisation by more than a unit in log evidence. Overall the highest 
evidences, over all datasets, were obtained with a power law spectrum containing a cutoff with 
a significant log evidence difference of roughly 2 units. The natural tilt and exponential cutoff 
present in the L-D spectrum is found to be favoured decisively by a log evidence difference 
of over 5 units, but only for a limited study within the best-fit concordance cosmology. 

Key words: cosmological parameters - cosmology:observations - cosmology:theory - cos- 
mic microwave background - large-scale structure 



1 INTRODUCTION 

The recent relea se of 3-year Wilkin son Microwave Anisotropy 
Probe (WMAP3: lffinshaw et al.|[2007h data have provided precise 
measurements of temperature fluctuations in the cosmic microwave 
background (CMB). The accepted inflationary paradigm suggests 
that a primordial spectrum of almost scale-invariant density fluctu- 
ations produced during inflation went on to produce the observed 
structure in the CMB and that seen on large-scales in the current 
distribution of matter. Now, for the first tim e a purely scale in- 
variant primordial spectrum is ruled out at la JSpergel et all2007l : 
IParkinson et alj |2006) in favour of a 'tilted' spectrum with n < 1. 
The WMAP team have al ready attempted l imited constraints on the 
form of the spectrum and lParkinson et al. ( 2006) have conducted a 
model selection study to ascertain the necessity of a tilt in the spec- 
trum with the new data. In this paper we extend both studies to a 
suite of models covering a wide variety of possibilities based on 
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both physical and observational grounds. We use a fully Bayesian 
approach to determine the model parameters and comparative evi- 
dence to ascertain which model the data actually prefers. 

Our previous paper teridges etaT]|2 006) [Bri dges06] used 
WMAP 1-year data (WMAP1; Ben nett et alj 120031) to constrain 
the same set of models. These generalisations were motivated 
principally by observations of a decrement in power on large- 
scales from WMAP1 and a tilting spectrum on small-scales from 
high resolution experiment s such as the A rcminute Cosmology 
Bolometer Array (ACBAR; iKuo et alj|2004) . the Very Small Ar- 
ray (VSA ; Dickinson et al. 2004) and the Cosmic Background Im- 
ager (CBI; Rea dhead et alj|2004h . With two more years observing 
time and improved treatment of systematic errors the decrement 
in power on large scales is now somewhat reduced, yet still evi- 
dent in WMAP3 and is now constrained almost to the cosmic vari- 
ance limit while the tilting spectrum on small-scales is now seen 
even without the aid of high-resolution small scale experiments, 
due to tighter constraints on the second acoustic peak. On phys- 
ical grounds we test a broken spectrum caused perhaps by dou- 
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ble field inflation dBarriga et al feoOll) and a spectrum predicted by 
lLasenbv & Doranl 1 20051) (L-D) naturally incorporating an expo- 



nential cutoff in power on large scales by considering the evolu- 
tion of closed universes out of a big bang singularity, with a novel 
boundary condition that restricts the total conformal time available 
in the universe. We also aim to reconstruct the spectrum in a num- 
ber of bins in wavenumber k. 



2 MODEL SELECTION FRAMEWORK 

Bayesian model selection is now well established within the 
community as a reliable means of appropriately determin- 
ing the most efficient parameteri sation for a model, penalis- 
ing any unecessary complication ijaffe 19961, Drell et al. 2000L 
John & Nar likar 2002, [Hobson et al j|200i iHobson & McLachlanl 
200l ISlosar et alj 120031 ISaini et alj l2004l iMarshall et al] 1 20031. 
Niarchou et alJ |2004 iBasset et all |2004 iMukheriee et alj 12000. 
Trottj 120071 beltran et alj I2005L Bridges06). Recently much 
progress has b een made in improving the speed and accuracy of ev- 
id ence results jParkinson et al. 2006) by implementing the method 
of lSkillingl d2004) known as nested sampling. In this paper we wil l 
employ our own implementation of this method (Sha w et alj2007t) 
to evaluate the evidence and use standard Markov Chain Monte 
Carlo (MCMC) to make parameter constraints. 



2.1 Markov Chain Monte Carlo sampling 

A Bayesian analysis provides a coherent approach to estimating 
the values of the parameters, ©, and their errors and a method for 
determining which model, M, best describes the data, D. Bayes 
theorem states that 



P(0|D,M) 



P(D|0,M)P(0|M) 
P(D\M) ' 



(1) 



where P(0|D,M) is the posterior, P(D[0,M) the likelihood, 
P(&\M) the prior, and P(D\M) the Bayesian evidence. Conven- 
tionally, the result of a Bayesian parameter estimation is the poste- 
rior probability distribution given by the product of the likelihood 
and prior. In addition however, the posterior distribution may be 
used to evaluate the Bayesian evidence for the model under consid- 
eration. 

We will employ a MCMC sampling procedure to explore the 
posterior distribution using an adapted version of the COSMOMC 
package dLewis & Bridlel2002h with four CMB datasets; WMAP3, 
ACBAR the VSA and CBI We al so include the 2dF Galaxy Red- 
shift Survey dPercival et alj|200ll) . the Sloan Digital Sky Survey 
( Abazajian et al. 2003 j) and the H ubble Space Telescope (HST) 
key project i Fre edman et al.l200ll) . These set of experiments com- 
prise dataset I. Additionally we include two datasets which cover 
different scales and probe independent source s. The Ly-a forest 
dMcDonald etal] 120061 : iMcDonald et alj|2005l) (which combined 
with dataset I makes up dataset II) comprises cosmological absorp- 
tion by neutral hydrogen observed in quasar spectra in the inter- 
galactic medium. It probes fluctuation scales that are small (~ Mpc) 
in comparison to the other datasets used at redshifts between 2-4 
so that primordial information has not been erased by non-linear 
evolution. It thus provides a very useful complementary observa- 
tion when c onstraining the form of the primordial spectrum. Previ- 
ous authors dViel. Haehnelt & Lewisl2006l : ISeliak et alj200^) have 
already examined this dataset in conjunction with WMAP3 and 
others and found that most of the interesting results observed by 



ISpergel et ail d2007h . namely lowered scalar amplitude and a non- 
vanishing running index can be removed wh en Ly-q is included. 
Observations of luminous red galaxies (LRG) iTegmark et al1l2006l 
(when combined with dataset I becomes dataset III) consists of > 
46,000 galaxies taken from the full SDSS catalogue which rep- 
resent a highly uniform galaxy sample containing only luminous 
early-types over the entire redshift ran ge studied cons t ituting an 
excellent tracer of large scale structure. ITegmark et alj 120061) re- 
cently detected baryonic acoustic oscillations in the matter power 
spectrum extracted from this dataset, providing a welcome confir- 
mation of early universe physics in large scale structure data. 

In addition to the primordial spectrum parameters, we param- 
eterise each model using the following five cosmological parame- 
ters; the physical baryon density Qth 2 ; the physical cold dark mat- 
ter density Q. c h 2 ; the curvature density f^; the Hubble parameter 
h (Ho = h x lOOkms" 1 ) and the redshift of re-ionisation z re . 



2.2 Bayesian evidence and nested sampling 

The Bayesian evidence is the average likelihood over the entire 
prior parameter space of the model: 



£(e c , e B )P(ec)P(e B )d Jv ecd M eB, 



(2) 



where N and M are the number of cosmological and Bianchi pa- 
rameters respectively. Those models having large areas of prior pa- 
rameter space with high likelihoods will produce high evidence 
values and vice versa. This effectively penalises models with 
excessively large parameter spaces, thus naturally incorporating 
Ockam's razor. 

The method of nested sampling is capable of much higher ac- 
curacy t han previous metho ds such as thermodynamic integration 
(see e.g. lBeltran et alfc OQS, Bridges06). due to a computationally 
more effecient mapping of the integral in Eqn.[2]to a single dimen- 
sion by a suitable re-parameterisation in terms of the prior mass X. 
This mass can be divided into elements dX = ^(0)^0 which 
can be combined in any order to give say 



X(\) 



n{&)d N &, 



(3) 



C{&)>\ 



the prior mass covering all likelihoods above the iso-likelihood 
curve C — A. We also require the function C{X) to be a singular 
decreasing function (which is trivially satisfied for most posteriors) 
so that using sampled points we can estimate the evidence via the 
integral: 



Z = 



C{X)dX. 



(4) 



Via this method we can obtain evidences with an accuracy 10 times 
higher than previous methods for the same number of likelihood 
evaluations. 

A standard scenario in Bayesian model selection would re- 
quire the computation of evidences for two models A and B. The 
difference of log-evidences \xiZa — In Zb, also called the Bayes 
factor then quant ifies how wel l A may fit the data when compared 
with model B. I Jeffreys! fl 961) provides a scale on which we can 
make qualitative conclusions based on this difference: AlnZ < 1 
is not significant, 1 < AlnZ < 2.5 significant, 2.5 < AlmZ < 5 
strong and Aln.Z > 5 decisive. 
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3 PRIMORDIAL POWER SPECTRUM 
PARAMETERISATION 

3.1 H-Z, power-law and running spectra 

The early Universe as observed in the CMB is highly homogeneous 
on large scales suggesting that any primordial spectrum of density 
fluctuations should be close to scale invariant. The H-Z spectrum is 
described by an amplitude A for which we assume a uniform prior 
of [15, 55] x 10" 8 . Slow-roll inflation, given an exponential poten- 
tial, predicts a slightly 'tilted' power-law spectrum parameterised 



P(k) = A 



(5) 



where the spectral index should be close to unity; we assume a 
uniform prior on n of [0.5, 1.5]; fco is the pivot scale (set to 0.05 
Mpc -1 ) of which n and the amplitude A are functions. For a 
generic inflationary potential we should also account for any scale 
dependence of n(k) called running, so that to first order: 



n-l+(l/2) ln(fc/fc )(dn/dlnfe) 



(6) 



where dn/d\nk is the running parameter n run ; for which we as- 
sume a uniform prior of [—0.15, 0.15]. 

The inclusion of WMAP3 in dataset I now places impressively 
tight constraints on all three models (see Fig.Q}. The unexpected re- 
duction in the value of the optical depth to reionisation, to r ~ 0.09 
has had the effect of reducing the overall amplitude of the power 
spectrum, due to the well known r-A degeneracy. This effect is no- 
ticeable in all cases but most particularly so in the H-Z spectrum 
in Fig.[TJ(c). For the first time the single spectral index model now 
exhibits a constraint, to la, of n„ = 0.95 ± 0.02 (see Fig. Q](b)) 
excluding the possibility of a scale-invariant spectrum at this con- 
fidence level. Furthermore, a pure power-law (with n run = 0) is 
also excluded at the la level with n run — —0.038 ± 0.030 (Fig.Q] 
(a)). 



The addition of Ly-a data further increases constraints on all spec- 
tral parameters (see Fig. |2j, particularly so on spectral running. 
While a scale invariant spectrum is also ruled out with this dataset 
(n s = 0.96±0.02) a running spectrum is not preferred, with a con- 
straint on n run — 0.015 ± 0.015 representing a doubling in accu- 
racy over dataset I. A further tension exists between the amplitude 
of fluctuations as found using dataset I alone and when combined 
with Ly-a, the latter preferring a mu ch larger value of a s and thus 
higher scalar fluctuation amplitude. Seliaketal. (2006) estimates 
this deviation to be at the 2a level and treat it as a normal statistical 
fluctuation and not a sign of some unaccounted systematic flaw in 
either dataset. Since two independent analyses of WMAP3 + Ly-a 
and WMAP3 + LRG both suggest no significant running, a con- 
servative conclusion is that the running observed with dataset I is 
simply a statistical anomaly albeit at close to 2a. 
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Figure 1. Marginalised parameter constraints with dataset I (black) for 
the H-Z, single-index and running models compared with our WMAP1 
(Bridges06) analysis (red). 
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Figure 2. Comparison of parameter constraints using dataset I (black) for 
the H-Z, single-index and running models with dataset II (red) and dataset 
III (blue). 



3.2 Large scale cutoff 

Confirmation by WMAP3 of the large-scale decrement in power 
at I ~ 2, close to the cosmic variance limit, supports the possi- 
bility of a spectrum with some form of cutoff. We reexamine the 
case of a sharp cutoff as did the WMAP team, parameterising the 



scale at which the power drops to zero, k c with a choice of prior 
[0.0,0.0006] Mpc -1 which was made to limit the study to regions 
up to i ~ 6 at which point an appreciable cutoff is no longer ob- 
served. 



P(k) 



0, 



(7) 
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Although this model P(k) is not continuous, as the physical spec- 
trum should be, it does give an upper limit on the average cutoff 
scale, which is useful when comparing for instance the L-D spec- 
trum which does predict a form for the cutoff. 

On small scales this spectrum behaves just as the single-index 
power law and so constraints on A and n are similar (see Fig. 
f3](. Our WMAP1 constraint on k c showed a non-zero likelihood 
for a cutoff at k = i.e. no cutoff. This likelihood is marginally 
lower with WMAP3 presumably due to the higher amplitude of the 
£ = 3 multipole which was more heavily suppressed in WMAP1 
observations. The peak in the likelihood encouragingly remains at 
around k c — 2.8 x 10 -4 Mpc -1 . The scale of this cutoff is much 
larger than anything probed by either Ly-a (in dataset II) or LRG 
(in dataset III), so neither dataset improves on this scale constraint 
(see Fig.[4]l. 



3.3 Broken Spectrum 

Multiple field inflation would produce a symmetry breaking phase 
transition in the early universe causing the mass of the inflaton field 
to change su ddenly, momentarily violating the slow-roll conditions 
jAdamsetal.1l 19971) . The resultant primordial spectrum would be 
roughly scale invariant initially, followed by a sudden break lasting 
roughly 1 e-fold before returning to scale invariance. Because the 
slow-roll conditions are violated it is not trivial to calculate the form 
of the break, however a robust expectation is that it will be sharp as 
the field undergo ing the phase tran sition evolves exponentially fast 
to its minimum tearriga et alj200ll) . We parameterise the spectrum 
as: 

(A, k^k s 
P(k) = <^ Ck a -\ k s <k^k e , (8) 
{ B, k > k e 

where the values of C and a are chosen to ensure continuity. Four 
power spectrum parameters were varied in this model: the ratio of 
amplitudes before and after the break A/B with prior [0.3, 7.2]; 
k a indicating the start of the break with prior [0.01, 0.1] Mpc -1 ; 
\n(k e /k a ) to constrain the length of the break with prior [0, 4] and 



Figure 5. Marginalised ID and 2D probability constraints for the broken 
spectrum model, for k s , ln(k e /k s ) and A/B for our WMAP1 (Bridges06) 
analysis (red) and dataset I (black). 2D constraints plotted with 1<t and 2<r 
confidence contours. 



normalisation A with prior [14.9, 54.6] x 10 -8 . These represent a 
very conservative set of priors that allow the model a large degree 
of freedom in both the position of the break (which could occur 
anywhere from k ~ 0.01, well above any possible large scale cut- 
off) and the form, which could be extended, so as to mimic a tilted 
spectrum or occur as a sharp drop. In addition a prior that k e could 
not exceed 0.1 Mpc -1 was imposed so that only a region well cov- 
ered by the datasets used was explored. 

In this parameterisation a scale invariant spectrum would have 
an amplitude ratio A/B = 1 and a large value of \n(k e /k s ). Nei- 
ther our WMAP1 analysis nor dataset I suggest this to be a plau- 
sible explanation (see Fig. [5} with a distinct drop in amplitude 
(A/B ~ 1.2) starting on scales below k s ~ 0.01 Mpc -1 . The 
bimodal distribution in \n(k e /k s ) vs. k s observed with WMAP1 
is preserved, though less pronounced, with dataset I. This implies 
a preference for both a sharp break at a single scale of k « 0.04 
Mpc -1 and an extended break begining at k ~ 0.01. A drop in 
power is clearly a crude way of approximating a tilted spectrum, 
so the extended break was an expected result. The sharp break, 
could be in dicative of a some early universe physics or it could 
simply be as lBridleetalJ<2003h suggest in their WMAP1 analysis, 
an artifact of the inclusion of the large scale structure datasets. The 
transition scale lies close to the point at which large scale structure 
data becomes statistically significant. Below this scale the WMAP 
data would disfavour a break and above it large scale structure data 
would. The inclusion of both Ly-a and LRG (see Fig.[6j data shifts 
the start of the break to larger scales (most pronounced with LRG), 
but simultaneously lowers the amplitude of the break, thus produc- 
ing a more gradual extended slope, providing a better approxima- 
tion to a tilted spectrum. Thus, our results do not conclusively sug- 
gest that a break does exist in the primordial spectrum, rather they 
seem only to confirm the need for a spectral tilt. 
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Figure 7. Reconstruction of the primordial power spectrum in 8 bands of k 
(with la errors) for dataset I (red) and our WMAP1 (Bridges06) analysis 
(yellow) compared with the best fit 1-year (dotted) and 3-year (filled) H-Z 
spectrum. 



In(k7k ) 



Figure 6. Marginalised ID and 2D probability constraints for the broken 
spectrum model, for k B , ln(k e /k B ) and A/B for dataset I (black), dataset 
II (red) and dataset III (blue). 2D constraints plotted with la and 2a confi- 
dence contours. 



3.4 Power spectrum reconstruction 

Many previous attempts have been made to reconstruct the 
primo r dial spectrum directly f r om da ta: Wang 1994; Bridl e et all 



2006; 
2004 ; 
2007) 



20031; IShafieloo & Souradeed |2004| ; ISinha & Souradeed 
Tocchini-Valentini. Douspis & Silkl 120051 ; iHannestad 
Bridges06. Most recently the WMAP team dSpergel et alj 
attempted to reconstruct it using a set of amplitude bins in k using 
WMAP3 data alone. The method suffers from the natural side 
effect of imposing correlations between neighbouring bins and 
broadening other constraints. Therefore searching for features by 
this method is unreliable and difficult. However from a model 
selection viewpoint it does allow full freedom for the data to 
decide how many parameters are required of the model. Our 
parameterisation linearly interpolates between eight amplitude 
bins a n in k on large scales between 0.0001 and 0.11 Mpc -1 
parameterised logarithmically with fc^+i = 2.75fci so that 



P(k) 



ki < k < k i+ i 
k <5 k n 



(9) 



As expected an obvious tilt is discernible in our dataset I recon- 
struction (Fig. [7} the mean bin amplitudes deviating significantly 
from the best fit scale invariant spectrum at high k, not observable 
in our WMAP1 analysis. Though, within la limits a H-Z spectrum 
can still be fitted to both sets of data, however it would require a 
lower amplitude to accommodate the values at large k. No cutoff is 
observed as was hinted at in our WMAP1 analysis; however uncer- 
tainty at this scale is dominated by high cosmic-variance, making 
constraints inherently difficult. We find the addition of LRG data 
in dataset III provides little further constraint beyond that obtained 
from dataset I alone (see Fig.[8jl. The effect of including Ly-a data 
in dataset II is however, marked. Firstly the overall amplitude is 
lifted, owing to the larger as required by Ly-a. Furthermore lit- 
tle obvious tilt is discernible in full agreement with our analysis in 
Sec. 13. II Encouragingly similar features are seen in all three spectra 
such as the peaks at 0.003 Mpc -1 and 0.05 Mpc -1 . 
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Figure 8. Reconstruction of the primordial power spectrum in 8 bands of 
k (with la errors) for dataset I (red), dataset II (magenta) and dataset III 
(green) compared with the best fit (black) H-Z spectrum from dataset I. 



3.5 Closed universe Inflation 



lLasenbv & Dorad d2005h arrived at a novel model spectrum by con- 
sidering a boundary condition that restricts the total conformal time 
available in the Universe, and requires a closed geometry. The re- 
sultant predicted perturbation spectrum encouragingly contains an 
exponential cuto ff (as previously suggested phenomenologically by 
Efsta thiod[2003h at low k that yields a corresponding deficit in 
power in the CMB power spectrum. The shape of the derived spec- 
trum for a cosmology defined by Slo = 1.04, flt,h 2 — 0.0224, h = 
: 0.110 was parameterised by the function: 



0.6, flcdmh : 

P(k) = A(l 
where y = In 



0.023y) 2 (l - exp(-(y + 0.93)/0.47)) 2 , (10) 



x 3 x 10'M > -0.93. However we do 



Ho/100 

know this form to be fairly stable to changes in cosmology, partic- 
ularly in the position of the cutoff. We therefore analyse this spec- 
trum using dataset I, varying only the amplitude within the best fit 
'concordance' cosmology found for the single index model in Sec. 
13.11 As expected the amplitude is reduced when compared with 
our WMAP1 analysis (A = 27.1 ± 0.2 Mpc -1 ) due to the revised 
WMAP3 value for r. At small scales the spectrum then lies roughly 
in line with a tilted spectrum with n 3 ~ 0.96 (see Fig.[9](. For com- 
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Figure 9. L-D spectrum (solid-red) shown with the WMAP3 best fitting H- 
Z (solid-black), single-index with a cutoff (dotted-orange), without (dotted- 
blue), a running index (dashed-pink) and broken (dotted-black). 



scale. As expected from the parameter constraints a running spec- 
trum is preferred, but only at a significant level if narrow priors 
are chosen on n s and n run . Large scale power is suppressed in 
WMAP3 as it was in WMAP1, with little improvement in uncer- 
tainty: both datasets are almost cosmic-variance limited on these 
scales so a cutoff is preferred by a significant margin of over 2 
units in log-evidence. The evidence in favour of a broken spectrum 
is likely only due to its mimicking of a tilted spectrum as opposed 
to modelling of an abrupt break. 

Inclusion of Ly-a in dataset II confirms the conclusions drawn 
in Section [3~Tl showing a significantly reduced evidence in favour 
of spectral running (with narrow priors), from a +1.2 with dataset I 
to -1.3, a significant evidence difference between the two datasets. 
Unsurprisingly all three data combinations exhibit significant evi- 
dences in favour of a spectral cutoff, due presumably in each case 
to the decrement in WMAP3 data alone. 



Table 1. Priors placed on the principal primordial spectral parameters. 



A s (xl0" 8 ) 


[15,55] 


n s (wide priors) 


[0.5,1.5] 


n s (narrow priors) 


[0.8,1.2] 


n r un (wide priors) 


[-0.15,0.15] 


n rU n (narrow priors) 


[-0.1,0.05] 



parison the plot also shows the other best fit spectra (including the 
bimodality in the broken spectrum). 



3.6 Model Selection 

We shall now turn to the fundamental inference. Which of the 
models considered best describes the current data? As with our 
WMAP1 analysis we perform the analysis in two stages to accom- 
modate the L-D spectrum: the first, within the full parameter space 
including all models bar L-D, while the second was carried out 
within the fixed, best-fit single-index cosmology (as determined us- 
ing dataset I) including L-D. While this artificial division may seem 
flawed it should be remembered that within this particular fixed 
cosmology the only spectrum that is favoured is that of the single- 
index model itself, any other models should only be penalised. We 
chose two sets of priors (see Table [TJ on n s and n run , the wider 
of which provides a very conservative range given the constraints 
available from current data. The results have been normalised to the 
single-index power law spectrum with wide priors. 



3.6.1 Full cosmological parameter space exploration 

Our corresponding analysis of WMAP1 data was unable to make 
any conclusive model selection; with statistical uncertainty domi- 
nating results. With WMAP3 data and the method of nested sam- 
pling however we have arrived at a cusp in cosmological model 
selection. With dataset I, Table [2] confirms the disfavouring of a 
scale-invariant spectrum that we saw in the parameter constraints 
in Sec. 13. II though not at a decisive level according to the Jeffrey's 



3.6.2 Primordial parameter space exploration 

The large preference in favour of a spectral cutoff and increasingly 
tight constraints on the universal geometry being marginally closed 
(a fact that is p a rticula rly reinforced when examining the LRG data 
iTegmark et al.l d2006f) ). suggests that the L-D spectrum may pro- 
vide a very good fit to current data. For the remainder of this anal- 
ysis we will examine the L-D spectrum using dataset I only. 

As with the WMAP1 analysis, dataset I, also prefers the L- 
D spectrum, now by a slightly larger log-evidence (see Table [3]), 
which according to Jeffreys' scale now constitutes a decisive model 
selection. In Bridges06 we concluded that the 'significant' model 
detection was due to the form of the spectrum on large scales i.e. 
its exponential cutoff. However on small scales it behaves much 
like a tilted spectrum with slight running, a form we now know fits 
the data very well. Could it be that these small scale features are 
driving this model selection? 

To test this hypothesis we have analysed three hybrid spectra, 
where we have divided the L-D spectrum about k ~ 0.008 Mpc -1 
(denoted by the arrow in Fig. [9) which corresponds loosely with the 
angular scale (Z = 8) where a cutoff ceases to be observed. We have 
then spliced both sections, about this point, with various single- 
index spectra at large or small scales. Model A combines a single- 
index spectrum below k = 0.0008 Mpc -1 with L-D thereafter; 
Model B: an exponential cutoff from L-D with a fixed single-index 
model (n = 0.94) for k > 0.0008 and Model C: as for B but with 
a varying single-index model. 

The results (see Table[3j bear out our assertion; models A and 
B have roughly identical log-evidence values to the original L-D, 
demonstrating the data to be essentially indifferent to the presence 
of a cutoff. In comparison, model C is typically just as good as a 
power law spectrum (i.e. an evidence close to 0). This suggests that 
the L-D spectrum is attractive to the data as it naturally incorporates 
a tilt without the need to parameterise it. However the tilt present 
in the L-D spectrum (n s ~ 0.96) does coincide well with the best 
fit values obtained in Sec. l3.1l for a power law spectrum (which are 
fairly invariant among the datasets I, II or III). This fact coupled 
with a significant evidence preference for a large scale spectral cut- 
off suggests the L-D spectrum does provide a uniquely good fit to 
current data. 
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Table 2. Differences of log-evidences with respect to single-index model (with the widest priors) using the full cosmological parameter space. [Note evidence 
comparisons can only be made between models using the same dataset.] 



Model 


datset I 


+ Ly-o 


+ LRG 


H-Z 


-0.7 ±0.3 


-0.1±0.3 


-0.4 ±0.3 


Single-Index (wide priors) 


+0.0 ± 0.3 


+0.0± 0.3 


+0.0 ± 0.3 


Single-Index (narrow priors) 


+1.0 ±0.3 


+0.2± 0.3 


+0.7 ± 0.3 


Running (wide priors) 


-2.9 ±0.3 


-1.6±0.3 


-1.8 ±0.3 


Running (narrow priors on n s ) 


+0.4 ± 0.3 


-0.7±0.3 


+1.7 ±0.3 


Running (narrow priors on n run &n B ) 


+1.2 ±0.3 


-1.3± 0.3 


+1.0 ±0.3 


Cutoff 


+2.3 ± 0.3 


+1.7± 0.3 


+2.9 ± 0.3 


Barriga 


+1.0 ±0.3 


+1.2± 0.3 


+0.9 ± 0.3 



Table 3. Differences of log-evidences (primordial parameters only) for all 
models with respect to the single-index models preferred cosmology: Ho = 
1.035, n b h 2 = 0.0221, h = 0.58,U cdm h 2 = 0.112, with our previous 
WMAP1 results for comparison. 



Model WMAP1 (Bridges06) dataset I 

Single-Index +0.0 ± 0.5 +0.0 ± 0.2 

L-D +4.1 ±0.5 +5.2 ±0.2 

L-D (A) - +5.0 ± 0.2 

L-D (B) - +5.2 ± 0.2 

L-D (C) - +0.9 ± 0.2 



4 CONCLUSIONS 

A scale-invariant spectrum is now largely disfavoured by the 
dataset I with a spectral index n 3 — 0.95 ± 0.02 deviating by 
at least 2a from n 3 — 0. Moreover a running spectrum (n run = 
—0.038 ± 0.030) is now significantly preferred but only using the 
most constraining prior. The addition of Ly-a forest data improves 
all constraints but does not alter the preferred spectral tilt greatly. It 
does however, along with LRG data, suggest a significantly smaller 
running index (n run = —0.015 ± 0.015, n r un = 0.01 ± 0.05) . 
This tension has previously been analysed by ISeliaket ail 12006) 
who conclude that such discrepancies, even though at the 2<j level 
are consistent with normal statistical fluctuations between datasets. 
A power law spectrum with a cutoff provides the best evidence fit 
in our full parameter space study with a significant evidence ratio 
of roughly 2 units across all three datasets. The similarity of this 
cutoff model with the L-D spectrum suggests the latter should also 
provide a very good fit. This is indeed borne out with decisively 
large evidence ratios within our limited primordial-only analysis. 
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